Optimization of ship hull forms by changing CM and CB coefficients to obtain optimal seakeeping performance

Ship design involves optimizing the hull in order to enhance safety, economic efficiency, and technical efficiency. Despite the long-term research on this problem and a number of significant conclusions, some of its content still needs to be improved. In this study, block and midship coefficients are incorporated to optimize the ship’s hull. The considered ship was a patrol vessel. The seakeeping analysis was performed employing strip theory. The hull form was generated using a fuzzy model. Though the body lines generated by the midship coefficient (CM) and block coefficient (CB) varied indecently, the other geometric parameters remained the same. Multi-objective optimization was used to optimize CB and CM. According to the results of this study, these coefficients have a significant impact on the pitch motion of the patrol vessel as well as the motion sickness index. Heave and roll motions, as well as the added resistance, were not significantly influenced by the coefficients of CM and CB. However, increasing the hull form parameters increases the maximum Response Amplitude Operator (RAO) of heave and roll motions. The frequency of occurrence of the maximum roll RAO was in direct relation with CB and CM. These coefficients, however, had no meaningful impact on the occurrence frequency of other motion indices. In the end, the CB and CM coefficients were selected based on the vessel’s seakeeping performance. These findings might be used by shipbuilders to construct the vessel with more efficient seakeeping performance.


Introduction
Seakeeping refers to a vessel's ability to navigate safely during long storms and withstand rough conditions at sea.Seakeeping performance can be described as the dynamic behavior of a vessel in wind, waves, and currents.Among the performance indicators, there are comfort, crew workability, and damage to the vessel and cargo resulting from slamming and green water, as well as the vessel's earning potential.
There are six degrees of freedom on each vessel, including three linear motions (surge, sway, heave) and three angular motions (roll, pitch and yaw).In the first steps of designing new vessels, the determination of these operational motions, either analytically or non-analytically, can lead to the development of design conditions.Seakeeping performance is dependent upon the shape of the hull, which is designed by naval architects to achieve the highest level of performance.The sea-going ship that operates in open waters rarely sails in calm weather.On the contrary, ship's behavior at sea is often affected by waves and wind.When a ship navigates in a seaway, the ship 0 s forward speed decreases, compared to that in calm sea, because of added resistance due to winds, waves, rudder angle, and so forth.The magnitude of added resistance is about 15-30% of calm-water resistance [1].The immediate effect of waves is ship's motions and accompanying phenomena, such as accelerations.Ship accelerations, in turn, particularly vertical ones, impact on the human body and may cause motion sickness.The term "motion sickness", on ships known as sea sickness, is understood as a sickness due to ship motions that results in physical discomfort, with such symptoms as irregular breathing, nausea, vertigo, paleness and vomiting.In extreme cases a passenger or crew member has to be transferred to hospital.The actual reason for sea sickness is lack of conformity between different stimuli, eye signals and the labyrinth (inner ear), received by the human brain.People mainly suffer from sea sickness under deck, where the eye does not register any stimuli that the labyrinth would interpret as motion [2].
Due to the abovementioned issued, the optimization of the hull is one of the most important aspects of ship design in order to improve both safety and economic and technical efficiency for ships.Despite the fact that this problem has been studied for a long time and has achieved many significant results, some contents need to be improved in practical applications.Among these are the presentation and transformation of existing hulls, the development of an optimal mathematical model for a particular type of ship, or the development of a method for solving the objective function, etc.
Optimization of a vessel hull form includes a number of nontrivial issues, including selection of an appropriate function, choosing optimization scheme, geometrical representation of hull surface and choice of related design variables and constraints, selecting a practical and robust numerical tool for evaluating the objective function, and decision to perform optimization for a single point design or for multiple point design, e.g. for a single ship speed or for a range of speeds [3].Different numerical schemes have been used for studying the ship hull form optimization problem.Some of the most relevant studies are listed in Table 1.
Maxsurf is a naval architecture software that provides integrated tools for hull modeling, stability, motions and resistance prediction as well as structural modeling.This software is a well-known numerical tool for ship analysis [53][54][55][56][57] that also used in recent ship optimization studies [58][59][60].This software was also used in previous studies of the authors.Khosravi Babadi and Ghaseemi (2013) performed a number of numerical and experimental studies, including investigating the effect of variations in some geometrical parameters and hull form coefficients, including water plane coefficient (C wp ) and prismatic coefficient (C p ), on sea-
While hull form optimization has been studied extensively, the effects of changing main hull form coefficients such as C B and C M on motions have not been investigated.This paper fills that gap.
Although some recent studies have addressed the issue of ship hull optimization, more research is required to examine the effect of different hull form coefficients, including ship midship and block coefficients.The purpose of this study is to provide insight on how to optimize the ship hull taking into account these two coefficients.This study involves the following steps being performed: 1. Assessing the effect of the C B and C M parameters of the vessel on seakeeping indices (heave, pitch, and roll motions, added resistance, and motion sickness).
2. Optimizing the vessel's form based on seakeeping parameters.
The key novelties of this work are: (1) using block and midship coefficients as optimization parameters to improve seakeeping and ( 2) developing a fuzzy model to generate hull form variations while keeping other parameters constant.
Unlike previous hull form optimization studies that use complex geometrical parameters, we use only C B and C M as variables in a novel fuzzy hull generation model to provide useful insights.
In order to achieve the above-mentioned goals, a fuzzy model was developed in order to generate the hull form.The body lines generated by C M and C B in this model vary indecently, but do not affect the other geometric parameters (C P ,L (Length) and B (Breadth)).An index of seakeeping performance (SPI) is defined as an objective function expressing dynamic behavior.Using multi-objective optimization, some values of C B and C M were derived that optimize seakeeping.
It will be shown that by optimizing these coefficients, pitch and MSI will improve.On the other hand, the effect of these coefficients on the roll and heave motion as well as the added resistance is negligible.

The mathematical procedure
In this section, the mathematical procedure made to optimize the vessel hull form is presented.

Strip theory
The estimation of ship motions in the presence of regular waves, arbitrary heading, and constant forward speed of the ship, as well as the calculation of wave-induced horizontal and vertical shear forces, bending moments, and torsional moments, are based on "strip theory".According to the theory, the coefficients of the related ship in the ship motion equations in two dimensions are calculated, and then integrated throughout the ship length and transformed into the 3D global coefficients.
The basic idea of the strip theory is dividing the hull into several slices along the longitudinal direction as shown in Fig 1 .Under given loading conditions and speed conditions, for any combination of wave frequency and wave direction, the hydrodynamic coefficients, such as additional mass, additional damping, Froude-Krylov wave force, and diffraction wave force, were calculated on each slice by applying a unit amplitude regular wave to the hull.Finally, the force of each slice is integrated longitudinally to obtain the force of the entire hull.In a regular wave, a ship's movement may be broken down into two separate issues to solve.
1. Radiation issues: Only the ship's free swing motion is taken into account because there is no impact wave.This condition's hydrodynamic force is made up of words for increased mass force, damping force, and restoring force.
2. Diffraction problem: Only the impacts of regularly occurring incident waves on the hull are considered, assuming that the ship is stationary.Wave forces make up the hydrodynamic force at this moment.Incident wave force and diffraction wave force make up wave force.The latter is the wave force produced by the wave when it contacts the hull, whereas the former just takes into account the impact of the incident wave on the hull and ignores the impact of the hull's presence on the flow field.When the ship motion responses are linear and harmonic, the coupled six degree of freedom (6-DoF) equation of motion in the frequency domain is as follows: where M jn and A jn are the generalized mass and added mass matrices, respectively, B jn and C jn are the damping and restoring coefficients.The hydrodynamic coefficients representing the jth degree of freedom caused b mentioning the nth degree of freedom.F j is the exciting force and moment.Z n ; _ Z n andÂ Z n are the displacement, velocity and acceleration of the nth DOF, respectively.F j and w e are the forces (or moments) acted to the ship and wave encounter frequency.
The frequency domain transfer function of the hull motion is then calculated by substituting the hydrodynamic and wave forces into the equation for the six-degrees-of-freedom motion of the hull.The relative motion connection makes it clear that the motion response may be acquired at any point along the hull and that time differentiation can be used to get the appropriate speed and acceleration.Then, the hydrodynamic and wave forces are substituted into the hull 6-DoF motion equation to obtain the frequency domain transfer function of the hull motion.It is known from the relative motion relationship that the motion response at any position of the hull can be obtained, and the corresponding speed and acceleration can be obtained by time differentiation.For the motion calculation, the total hydrodynamic coefficients were computed with the Salvesen-Tuck-Faltinsen (STF) strip theory that is well known, so the details are not described here.
Response amplitude operator (RAO).RAO is also known as the transfer function (similar to the response curve of an electronic filter), describes how the response of the vessel changes with frequency.These are usually dimensionless due to the height of the waves and the slope of the waves.RAOs tend toward unity at low frequencies; this is where the ship simply moves up and down in waves and acts as a stopper.At high frequencies, the response tends to zero due to the effects of many destructive micro waves along the length of the ship.Usually, ships will also have larger peaks than unity; this occurs near the natural period of the circuit.The peak is due to resonance.An RAO value greater than unity indicates that the ship's response is greater than the amplitude (or slope) of the wave.
Motion response theory in irregular waves.We typically believe that the response of a ship's linear system may be superimposed homogeneously when predicting a ship's seakeeping performance in irregular waves.Additionally, the output is treated as a stationary random process when the input is a stationary random process.Under such assumptions, hydrodynamic calculation also known as the transfer function, may be used to determine the relevant connection between the response variable and the wave frequency (or period, wavelength) for each wave direction, each wave speed, and each loading condition (or response amplitude operators).The following formula may be used to get the response spectral density function from the transfer function and wave spectral density function: where H(w,β,U,C) is the response amplitude operator transfer function, S R is the response spectral density function, β,U,C represents the heading angle, ship speed, and ship loading condition, respectively.w; H S ; T Z are the wave frequency (rad/s), significant wave height (m), and average zero-crossing period (s), respectively.The encounter frequency will alter with the heading angles as the ship moves across the waves.The link between the encounter frequency and the wave frequency is as follows because the wave spectral density function at the encounter frequency and the wave spectral density function at the wave frequency have the same amount of energy [76]: The motion variance is given by the area under the motion energy spectrum as Hence, the ffi ffi ffi ffi ffi ffi ffi m o e p represents the RMS motion, and the significant motion amplitude is twice the RMS motion.In addition, the RMS velocity and acceleration are given by ffi ffi ffi ffi ffi ffi ffi m 2 e p and ffi ffi ffi ffi ffi ffi ffi m 4 e p [76].

Motion sickness index (MSI).
Despite scientific observations and research, no exact relations have been determined between ship motions and motion sickness.McCauley and O'Hanlon estimated quantitatively the impact of ship motions on the percentage of people that would suffer from sea sickness.It turned out that vertical accelerations in particular were responsible for motion sickness, while rolling and pitching had slight influence.Additionally, it was found that at a frequency of 0.167 Hz the occurrence of motion sickness increased.The Motion Sickness Incidence (MSI) index is commonly used for assessing possible occurrence of the illness: where MSI is the motion sickness incidence index, erf is the error function, a v is the mean value of vertical accelerations at a selected point and μ MSI = −0.819+2.32(log 10 w E ) 2 .Optimization procedure.The flowchart of the optimization procedure is given in Fig 2.

Modeling fuzzy changes in desired parameters
Fuzzy logic has advantages in ship hull optimization, including: • Handling imprecise and uncertain data: Fuzzy logic effectively deals with imprecise and uncertain information, common in ship design and optimization.
• Flexibility: Fuzzy logic incorporates subjective human knowledge and expertise into the optimization process, making it adaptable to different design requirements.
• Non-linear relationships: Fuzzy logic models complex, non-linear relationships between design parameters and performance criteria, which traditional optimization methods may struggle with.
• Robustness: Fuzzy logic-based optimization methods are often more stable, handling variations and uncertainties in design parameters without significantly impacting results.
This paper presents a fuzzy structure model for midship coefficient changes.The variation of these coefficients does not affect the other geometric parameters.Changes should be applied in such a way as not to affect the volume of the vessel.It is, therefore, necessary to develop a mathematical model of the changes that will cause the exact change in the desired parameter while maintaining other geometric parameters constant while holding constant the volume.
A fuzzy membership function, which is an extended Bell membership function, is defined as a function that makes changes to the body's lines.As a result, if any other changes are made, it increases and decreases the same volume, thus maintaining a constant volume change.Bell function is expressed as follows and its shape is a modified Gaussian distribution.Fig 3 illustrates the distribution shape of the Bell function.
Fig 3 illustrates how the fuzzy function at the upper level makes good changes to the body line in the middle and returns the variant body lines to the original lines at both ends with instant cuts in both directions, with the bell acting as a coefficient of variation.It is now necessary to restore the taken area of one side to the opposite direction using the same function used in the interval (0, -1).Therefore, this function is multiplied by a sine function.The fuzzy function coefficients in this model are defined as follows: The following is a description of the effects of parameters on the function: Coefficient b: Coefficient b is a positive number and is usually considered to be 4.The number corresponds to the upper part of the curve and contributes to the flat portion of it.This coefficient has a significant impact on the performance of the function.Coefficient a: Based on mathematical calculations and modeling, a coefficient equals 1.44 times the length of each line within a body.There will obviously be a difference in the coefficient of the body lines of the vessel.
Coefficient c: It determines the center of the bell function that is considered equal to the half-length of each line.
The sine function of the fuzzy function multiplied by Bell is given below: Data set x represents the length of the body lines to a point at which the fuzzy shift applies the same area to the other side.The length of the new points is determined by the function H (x). Using coefficient L at the end of the function, the different sizes of the changes are applied manually.This is the fuzzy function -neurological function obtained from Bell's model to apply to the vessel model: Using the obtained membership function, a fuzzy inference system is formed using a membership function H(x) and a simple inference function.Fuzzy relations are applied to any of the vessel body lines based on the customization done.As a result of using the above model, the vessel volume is not altered and the coefficient C B and C M are affected positively by adjusting geometrical values.

The objective functions in the frequency domain
Based on the original model, 10 variant models are generated for coefficients C B and C M and seakeeping parameters which include roll, pitch, heave, added resistance, and motion sickness acceleration calculated at three heading angles.For an irregular wave, the dynamic behavior of a vessel is calculated using the strip theory method.As the coefficient changes are very small, the results of calculations are very similar for different forms of the body.Some wave frequencies, however, demonstrate differences in dynamic behavior.Most of the difference can be found in the areas of maximum and minimum amplitudes.As a result, the objective of optimization is to reduce maximum amounts, thereby reducing the effective amounts of seakeeping parameters.Five objective functions, including roll, pitch, heave, seasickness acceleration, and added resistance in the waves, are plotted.As a result of the independent solution of the roll equation in the calculation, we will not have roll for head waves.For C B and C M , each diagram shows a variation range of ±3%.
In multi-objective optimization with seakeeping purposes, for each generated geometric model and every heading angle, five RAO curves will be generated based on the seasickness accelerations.These curves will include roll, pitch, heave, added resistance, and MSI.The optimization does not take into account a certain frequency; the significant value for every curve or parameter within a certain frequency range will be calculated.For example, C B and C M coefficients with percentage changes of +3%, +2.5%, 2%, 1.5%0%, -1.5%, -2%, -2.5% and -3% and heading angle of 150 degrees can be considered nine significant for each parameter.
For the purpose of achieving the objective functions, polynomial curve fitting is used, contributing these nine points.The objective function can be applied to each coefficient (variable) at each heading angle.Each heading angle should be combined with weights for any parameter (such as roll) to determine the objective function.

Wave spectrum and motion indices spectrum
The vessel is planned to perform operation in the Gulf of Oman.For this sea, the Pierson Moskowitz wave spectrum is suitable to be used as the wave spectrum.Generalized Pierson Moskowitz spectrum is as follows: In this study, we used ITTC spectrum.For ITTC spectrum the input parameters are H S and one of the T E ; T p ; � T or T Z .Also, A and B parameters are defined as follows: Pierson-Moskowitz table for sea states (Beaufort) 3 and 5 are listed in Table 2: RAO is the relation between the vessel motions to the wave amplitude which is usually draws in the shape of a dimensionless diagram based on the incidence frequency or dimensionless frequency.For linear and angular motions, RAO is defined by Eqs 16 and 17, ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where n is the desired degree of freedom and K is the wave number.

Mathematical modeling of the main body
This step involves entering the body line points into the MATLAB database.The training is carried out using a neural network with a transfer function in the hidden layer and a linear transfer function in the output layer, and 29 neurons in the hidden layer with a training function.
Levenberg-Marquardt algorithm.This Levenberg-Marquardt algorithm combines gradient reduction algorithms with Gauss-Newton algorithms (GNA).Unlike the Gauss-Newton algorithm, the Levenberg-Marquardt algorithm often finds a solution even if it begins far from the final minimum.Updating the parameters of this algorithm is performed by the following formula: and gradient: where: In order to update the weights, the following equation is used: Jacobian matrices are calculated using the same procedure as gradient matrices, except that derivatives are used instead of differences.The algorithm of this method, assuming matrix X is the input matrix, is as follows: 2. For each hidden layer and output layer, calculate the Jacobian matrix: 3. By choosing the μ value, the weight change is calculated: 4. Upon determining the amount of performance change, the error is checked, and if it increases, μ is decreased until the error decreases.
5. The weights have been updated and we have returned to the first stage.

Seakeeping and seakeeping performance index (SPI)
The science of seakeeping involves investigating and predicting vessel motions.Designers are required to provide information regarding the vessel's seakeeping performance, including local motions and accelerations, added resistance, deck wetness, and bow slamming."Vessels suitable for seagoing" must be capable of maintaining their motion in harsh conditions, so that deck wetness does not occur (Fig 4).The vessel should be independent of wind or wave direction, continue to follow its intended course even in inappropriate locations, and be able to quickly adjust to small angle changes.In addition, it must maintain a steady speed without slamming or abnormal fluctuations in the torque of the power transmission axis.Motions of the ship must be within an acceptable range and vibrations should not be excessive.The Seakeeping Performance Index (SPI) is a common measure of how well a ship handles rough seas.It calculates the percentage of time that the ship stays within specific motion limits.The SPI depends on assumptions about the frequency of different sea states and the likelihood of different ship speeds and headings.To evaluate the SPI, we predict how the ship will move for each combination of heading, speed, and sea condition.We then compare these predictions to a set of criteria that determine the optimal performance limits for the ship's mission in that particular sea condition.

Simulation procedure
For changes to be applied to one parameter and other coefficients to remain constant, changes should be applied to the area under the draft, and anti-changes should be applied to the same area of the body lines.This will lead to the rise and fall of the center of gravity and subsequently the change of KB.For this purpose, the focus of the modeling is placed on the front lines of the vessel, which are called "tearing" lines.The modeling steps are introduced by the following steps: 1. Using the cross sections of the vessel's front lines, a curved line is selected.

Body line Fr13 divides the desired float into the upper part of the chine and the lower chine,
each of which is regulated using a polynomial matching method.Based on the results, we have obtained the following relationships: f(x) = 2.167x−3642 above the chain f(x) = 0.0002008x 2 +0.3912x+63.19bottom of chain 1.As a result of the trial-and-error method, we obtain a calculated model that is equal to the displacement area.
2. Polynomial matching is performed separately for each part of the modeled graph.

Effective seakeeping parameters
This study investigates the dynamic effects of the vessel in the wave, as well as the operability criteria according to the type of operation and sub-systems required by the vessel.The roll, pitch, yaw, seasickness acceleration, and added resistance at the headings of 120, 150, and 180 are of the utmost importance.In addition to satisfying many criteria, by reducing the fluctuation range of these parameters, the derivatives (speed and acceleration) of these parameters are also reduced.These parameters will be used to determine the seafaring performance index.RAO and effective seakeeping parameters.Following the process described in the previous sections, geometric coefficients were determined, and based on the range considered (±3%), 64 models were produced.In this section, we will calculate the relationship between the effective seakeeping parameters and the geometric coefficients.This means that when considering the effective seakeeping parameters, such as roll, pitch, heave, seasickness acceleration and added resistance at different wave impact angles, the goal is to obtain these parameters according to geometric coefficients.The following steps will be followed: 1. Calculation of the RAO spectrum for each coefficient in all three incident angles.
There is a RAO spectrum diagram for each parameter (e.g., Heave motion), each incidence angle, and each coefficient (e.g., C B ).Note that for the parameters of roll, pitch, heave and added resistance, the amplitude of wave excitation is used to calculate RAO.There is no RAO associated with motion sickness acceleration.
2. RMS calculation for each RAO spectrum curve.
3. Calculating effective parameters in terms of coefficients at each angle of incidence by normalizing the values.
An example of C M and C B modifications.It is necessary to change the points of each of the lines in order to apply the desired changes to the vessel's lines.We apply the changes in the set of points in the desired direction and analyze the results as a new modeled vessel.In this case, the systematic changes are applied in accordance with the appropriate fuzzy function.In this problem, the constant weight of the vessel is assumed, and the change in the variable is the According to the variable relationship of the C M coefficient, the goal is to maintain the numerator while changing the denominator.Under the waterline of the main body, for example, an area of the lower half of the line should be tilted to the right, and an equivalent area of the upper half should be tilted to the left.These changes will result in an increase in the output of the fraction's denominator and a decrease in the coefficient.
An example of the output of the proposed method for C B modeling is shown in Fig 6 .This figure shows the main lines of the body in blue and a designed model in red.

Choosing the best answer with fuzzy-promethea decision-making approach
In order to select the optimal values, the fuzzy-promethea approach is used.This method is a fuzzy multi-criteria decision-making method that performs the evaluation process by using the ranking of weighing coefficients.The ranking of options is done by comparing the pairs of options in each index.The comparison is measured based on a pre-defined superiority function with range [0,1].The superiority function p, for comparing two options a and b in terms of index j, is defined as follows: P j ða; bÞ ¼ P j ½d j ða; bÞ� ð20Þ where d is the superiority function in the promethea method.Six most common superiority functions are shown in Fig 7.
The final ranking with the priority of the two options is obtained by summing the priority of all indicators, which is called the final ranking or overall value and is obtained by the following relationship: pða; bÞ ¼ X k j¼1 w j p j ða; bÞ; ð where w j is the weight of the j-th index.If the number of options (n) is more than two, the final ranking is obtained by summing the values of pairwise comparisons.This also specifies the ranking order (input / output / final): ;ðaÞ ¼ ; þ ðaÞ À ; À ðaÞ ð24Þ

An overview of the ship under consideration
A single-hulled patrol vessel with a V-shaped hull is being considered, which moves at a speed of 35 knots.The main specifications of the ship can be found in Tables 3 and 4.These tables may be amended as long as the specifications and body shape have not been finalized.A 3D view of the vessel is shown in Fig 8.

Validation
In this section, the results obtained using numerical model are compared with the laboratory data.A frame of the experimental study is shown in Fig 9 .The wave amplitude is 1.5 cm.In this case, model presents more extreme motions in higher sea states.When the wave length is almost the same of the vessel length, the biggest heave and pitch motions are observed.This is evident in Figs 10 and 11.As these figures show, there is a suitable correspondence between experimental and numerical data.Hence, the numerical model is accurate enough to be used in further numerical investigations.

The initial state of seakeeping of the vessel under consideration
This section presents the response spectrum at three headings of 120, 150, and 180 degrees at 15 knots for a patrol vessel.

Spectrum of effective seakeeping parameters for C M coefficient
In the single parameter C M spectrum, RAO diagrams were calculated for heave, pitch, roll motions, added resistance, and MSI against the incidence frequency at a speed of 35 knots for different headings.

Spectrum of effective seakeeping parameters for C B coefficient
The RAO diagrams for heave, pitch, roll, added resistance, and MSI were calculated for headings of 150 degrees in the single-parameter C B spectrum at a velocity of 35 knots.Also included    RAO increased with a decrease in C B , on the other hand.Furthermore, there was no significant impact of C B changes on the frequency of occurrence of maximum added resistance RAO.
A significant impact of C B on MSI can be seen in Fig 26 .It was found that C B changes had a direct correlation with the maximum MSI RAO.C B changes did not significantly affect the frequency of occurrence of the maximum added resistance RAO.

The final optimal solution of the geometric parameters of the vessel
In this section, a final optimal prioritized table is given for each parameter, which are numbered in order of superiority of the answer.In this way, the best answer of the first rank will be created, and this importance will decrease for the answers of lower ranks.Tables 5 and 6 show the final optimal values for midship and block coefficients, respectively.

Conclusions
This study is provided optimization of ship hull form taking into consideration the block and midship coefficients (C M and C B ).To generate the hull form, a fuzzy model was developed.In this model, there is an indecent variation in the body lines generated by C M and C B , but this does not affect the other geometric parameters.An index of seakeeping performance (SPI) measured the dynamic behavior of the vessel.C B and C M were optimized using multi-objective optimization.These are the most important findings of this study: • The results reveal C B and C M have a significant impact on pitch and MSI which has not been shown before.The pitch motion can be uncomfortable for passengers and crew and can also affect the stability and performance of the vessel.Hence, decreasing this motion is of great importance in ship design.
• Roll motion was very little affected by both C B and C M coefficients.The maximum roll RAO was increased by increasing these coefficients, but decreasing C B had no meaningful impact on it.Furthermore, a direct relationship was found between changes in these coefficients and the frequency of occurrence of the maximum roll RAO.
• Heave RAO was not significantly impacted by block and midship coefficients.It was found that changes in these coefficients directly impacted the maximum heave RAO.Additionally, changes in these coefficients had no significant impact on the frequency of maximum heave RAO.
• Both coefficients considered had a significant impact on pitch motion.The change in these coefficients was found to have an adverse relationship with the maximum pitch RAO.There was, however, no significant effect of changes in these coefficients on the frequency of occurrence of the maximum pitch RAO.
• The coefficients for the block and midship did not have a significant impact on added resistance.This coefficient was found to have an adverse relationship with the maximum added resistance RAO.In addition, these coefficients had no significant impact on the frequency of occurrence of the maximum added resistance RAO.• It was found that both coefficients had a significant impact on the MSI.A direct correlation was found between coefficient changes and the maximum MSI RAO.It was also observed that the frequency of occurrence of the maximum added resistance RAO was not significantly affected by coefficient changes.Decreasing the MSI improved the conformability of the crew onboard the vessel.
• Based on the mathematical procedure and performed optimization, the optimum values for block and midship coefficients were evaluated.

3 .
Fig 5 shows the modeled body lines.

Fig 10 .
Fig 10.Pitch RAO of the original model at Fr=0.21 in head waves.https://doi.org/10.1371/journal.pone.0302054.g010 Figs 17-21 illustrate the RAO diagrams for heave, pitch, and roll motions, as well as added resistance and MSI.Figs 22-24 illustrate the RMS diagrams for heave, pitch, and roll motions.The Fig 17 illustrates the heave RAO for changing CM coefficients at different wave frequencies.There was no significant impact of C M on heave RAO.However, the maximum heave RAO was increased by increasing the C M .In contrast, decreasing C M resulted in a decrease in the maximum heave RAO.Further, the frequency of occurrence of the maximum heave RAO was not significantly affected by changes in C M .

Fig 18
Fig 18 illustrates that C M has a very small effect on roll motion.Increasing C M , however, increases the maximum roll RAO.Conversely, decreasing C M did not have a meaningful impact on maximum roll RAO.Furthermore, an increase in C M resulted in a decrease in the frequency of occurrence of the maximum roll RAO, whereas a decrease in C M resulted in an increase in the frequency of occurrence of the maximum roll RAO.According to Fig 19, C M had a significant impact on pitch motion.By increasing C M , the maximum pitch RAO was decreased.In contrast, decreasing C M increased the maximum pitch RAO.Moreover, changes in C M had no significant effect on the frequency of occurrence of the maximum pitch RAO.In accordance with Fig 20, CM had no significant impact on added resistance.As C M was increased, the maximum added resistance RAO was decreased.Conversely, a decrease in C M increased the maximum added resistance RAO.Also, C M changes did not significantly impact the frequency of occurrence of the maximum added resistance RAO.

Fig 21
Fig 21 illustrates that C M had a significant impact on MSI.There was a direct correlation between C M changes and the maximum MSI RAO.Furthermore, the frequency of occurrence of the maximum added resistance RAO was not significantly affected by C M changes.

Fig 24
Fig 24 shows that C B significantly influenced pitch motion.The maximum pitch RAO was decreased by increasing C B .Conversely, a decrease in C B resulted in an increase in the maximum pitch RAO.Additionally, changes in C B did not have a significant effect on the frequency of occurrence of the maximum pitch RAO.According to Fig 25, C B had no significant impact on resistance added.Increasing C B resulted in a decrease in maximum added resistance RAO.The maximum added resistance